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Irreversible processes of one-dimensional quantum perfect Lorentz gas is studied on the basis of the 
fundamental laws of physics in terms of the complex spectral analysis associated with the resonance 
state of the Liouville-von Neumann operator. A limitation of the usual phenomenological Boltzmann 
equation is discussed from this dynamical point of view. For a Wigner distribution function that 
spreads over moderately small scale comparative to the mean-free-path, we found a shifting motion 
in space of the distribution that cannot be described by the hydrodynamic approximation of the 
kinetic equation. The mechanism of the shifting motion has two completely different origins: one is 
due to different value of the imaginary part of the eigenvalue of the Liouvillian and predominates in 
moderately small wavenumber associated to the spatial distribution, while the other is due to the 
existence of the real part of the eigenvalue associated to a wave propagation and predominates in 
moderately large wavenumber. 

PACS number (s): 05.30.-d, 05.60.-k, 05.20.Dd 

I. INTRODUCTION 

The purpose of this paper is to derive an irreversible 
kinetic equation through the fundamental laws of physics 
described by the Liouville-von Neumann equation with¬ 
out relying upon phenomenological arguments. The ex¬ 
ample we discuss in this paper is a weakly-coupled one¬ 
dimensional (ID) quantum perfect Lorentz gas. 

The world surrounding us is irreversible, and there is 
an arrow of time oriented toward our future that breaks 
time-symmetry. Nevertheless, the fundamental equations 
of motion in microscopic level are time-symmetric. To re¬ 
solve this apparent contradiction is a long-standing prob¬ 
lem in dynamics and statistical mechanics. A traditional 
assertion found in many textbooks in statistical mechan¬ 
ics is that the irreversibility is not a fundamental dy¬ 
namical property, and it comes from our approximation 
of smearing out of detailed microscopic information that 
cannot be controlled due to our limitations as human be¬ 
ings (see e.g. 0 )- This anthropocentric view to explain 
the irreversibility is called coarse-graining. However, af¬ 
ter we found the dissipative structure that allows our 
existence, this anthropocentric view is incomprehensible, 
since we are not fathers of the arrow of time, but we are 
children of the arrow of time 00 - 

With Ilya Prigogine, one of the authors (TP) in this pa¬ 
per has formulated the irreversible dynamics through the 
fundamental laws of physics without relying upon such 
anthropocentric phenomenological arguments 0, 0• In 
this formulation we consider the spectral representation 
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of the Liouville-von Neumann operator (Liouvillian) for 
systems in which irreversibility is expected such as the 
case of the thermodynamic limit where intensive vari¬ 
ables and extensive variables exist. This formalism is an 
extension of the well-known Brillouin-Wigner-Feshbach 
formalism of the eigenvalue problem of the Hamiltonian 
in quantum systems to the eigenvalue problem of the Li¬ 
ouvillian. In this formulation, we start with a complete 
set of projection operators (often we choose a complete 
set of projection operators consisting of the eigenstates 
of the unperturbed Liouvillian). Then, for each subspace 
associated to a projection operator we construct an effec¬ 
tive Liouvillian that shares eigenvalues with the original 
Liouvillian, similarly to the case for the effective Hamil¬ 
tonian in the Hamiltonian dynamics. 

Corresponding to the self-energy part of the effective 
Hamiltonian, we have self-frequency part of the effec¬ 
tive Liouvillian. By the same reason in the effective 
Hamiltonian, the effective Liouvillian may become a non- 
Hermitian operator for unstable systems with resonance 
singularity, in spite of the fact that the Liouvillian itself 
is a Hermitian operator indicating the time-symmetric 
nature of the fundamental laws of physics. As a result, 
an eigenvalue of the Liouvillian becomes a complex num¬ 
ber which is located at a resonance pole of the resolvent 
operator (or, the Green’s function) of the Liouvillian in 
the complex frequency space. Associated to a complex 
eigenvalue, we have a resonance state for the Liouvil¬ 
lian as in the Hamiltonian case. For the Hamiltonian 
case, the imaginary part of a complex eigenvalue gives 
a decay rate of an unstable excited state. Similarly, the 
imaginary part of a complex eigenvalue of the Liouvillian 
gives a transport coefficient associated to an irreversible 
process. 

It is important to recognize that the existence of the 
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resonance pole of the resolvent operator of the Hamilto¬ 
nian is a direct consequence of the mathematical struc¬ 
ture of the Hamiltonian, and does not result from an an¬ 
thropocentric operation, such as the coarse-graining. In 
other words, the decay rate that breaks time-symmetry in 
the Hamiltonian system is an intrinsic dynamical prop¬ 
erty of the system. For exactly the same reason, irre¬ 
versibility associated with transport coefficients in the 
Liouvillian system is an intrinsic dynamical property of 
the Liouvillian dynamics derived from the microscopic 
fundamental laws in physics. 

In spite of the parallelism between the Hamiltonian 
case and the Liouvillian case, there is a fundamental dif¬ 
ference between them in the self-energy or self-frequency 
part in the effective operators mentioned above. Indeed, 
even in the case where the unperturbed Hamiltonian does 
not degenerate, the unperturbed Liouvillian has intrinsic 
degeneracy. For quantum case this is a result of the fact 
that the Liouvillian is defined by the commutation rela¬ 
tion with the Hamiltonian (see (|8]l). As a result, even in 
the case the self-energy part of the Hamiltonian is just 
a number, the self-frequency part becomes an operator. 
In statistical mechanics this self-frequency part operator 
has been called the collision operator, and it is a central 
object in non-equilibrium statistical physics. 

Even though main framework of our formalism men¬ 
tioned above of the irreversible processes is settled on 
the microscopic laws of the fundamental dynamics, its 
application to actual systems is generally extremely com¬ 
plex, because the infinitely many degrees of freedom are 
mutually interacting in systems that show irreversibility. 
Moreover, the kinetic equation derived from the eigen¬ 
value problem of the Liouvillian has a nonlinear struc¬ 
ture that comes from two different origins; (i) one is due 
to the appearance of the product of the one-body dis¬ 
tribution function due to the many-body effect, and (ii) 
the other is due to the fact that the collision operator 
itself depends on its eigenvalue. Due to these complex¬ 
ity, we have applied this formulation so far only to the 
case where drastic approximations to solve the complex 
eigenvalue problem of the Liouvillian is applicable, such 
as the cases, for example, near equilibrium state where 
a linearized approximation with respect to the product 
of the distribution function is applicable, or the spatial 
inhomogeneity has a macroscopic scale at least of the 
order of the mean-free-length [6H8| . Since one cannot an¬ 
alyze the dynamics in a microscopic scale much smaller 
than the mean-free-length by this drastic approximation, 
we have not yet been able to demonstrate the advantage 
of our microscopic approach over the phenomenological 
coarse-graining approach. 

The main purpose of this paper (and successive forth¬ 
coming papers) is to analyze microscopic origin of irre¬ 
versibility that is applicable to all spatial scale of the 
inhomogeneity including microscopic scale of the inter¬ 
action range for the weekly-coupled ID quantum perfect 
Lorentz gas where the mass of a light test particle is neg¬ 
ligible as compared with mass of heavy particle. The 


advantage of the Lorentz gas is a mathematical simplic¬ 
ity. Thanks to the disparity in number of degrees of free¬ 
dom between a light test particle and background many 
heavy particles, one can neglect the time evolution on 
the distribution function of the heavy particles in the 
thermodynamic limit. As a result, the complexity of the 
nonlinearity associated to the origin (i) mentioned above 
does not exist in the Lorentz gas. 

We especially focus our attention in this paper on lim¬ 
itations of the usual phenomenological Boltzmann ap¬ 
proximation. Through our microscopic approach based 
on the complex spectral representation of the Liouvillian, 
the meaning of the coarse-graining will become clearer. 

In addition, we report our discovery on an interesting 
shifting motion in space that is generated by the Boltz¬ 
mann collision operator, especially for the distribution 
function localized in moderately small spatial scale less 
than the mean-free-length, but yet with a large enough 
width as compared with the microscopic scale given by 
the interaction range. We found that there are two 
completely different mechanisms of the shifting motion 
in space; one is due to the asymmetry of the momen¬ 
tum distribution and predominates in a moderately small 
wavenumber associated to the inhomogeneity, while the 
other is due to the wave propagation associated to the 
real part in the eigenvalue that appears only in a mod¬ 
erately large wavenumber. As we shall see, this shifting 
motion cannot be described in the hydrodynamic approx¬ 
imation (such as the diffusion equation) that is applicable 
only to situations with extremely larger spatial scale as 
compared with the mean-free-length. 

The structure of the present paper is as follows: In 
Sec.II, we introduce the ID quantum perfect Lorentz gas. 
In Sec.Ill, we summarize general formalism of the com¬ 
plex spectral representation of the Liouvillian and apply 
it to the model. In Sec.IV, we discuss limitations of the 
Boltzmann approximation. In Sec.V, we present a solu¬ 
tion of the eigenvalue problem of the Boltzmann collision 
operator and discuss its properties. In Sec.VI, we discuss 
how the spectrum of the Boltzmann collision operator 
affects transport processes by analyzing time evolution 
of the Wigner distribution function. In the last section 
VII, we give concluding remarks. In Appendix, we show 
correction terms to the collision operator of the perfect 
Lorentz gas in a case where the ratio of the mass of the 
light particle to that of the heavy particles is not negli¬ 
gible. 


II. SYSTEM 

We consider a weakly-coupled ID quantum Lorentz 
gas. The Lorentz gas consists of one light-mass particle 
(the test particle) with mass m and N heavy particles 
with mass M. The Hamiltonian of the system is given 
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by 


H = H 0 +gV 



N 2 N 

+ E^+sE F ( \ x ~ x i\), ( 1 ) 

3 =1 3=1 


where g is a dimensionless coupling constant and the in¬ 
teraction V is assumed to be a short-range repulsive po¬ 
tential. In this paper, we consider a weak-coupling situ¬ 
ation (5 C 1). In the following analysis, we restrict our 
attention to the limit 


in which the system is called the perfect Lorentz gas |H, 

E3- ' 

We assume that the system is enclosed in a large ID 
box of volume L with the periodic boundary condition. 
The interaction potential is expanded in the Fourier series 
as 


V(\x-x j \) = ±J2 V 0’> eign{X ~ Xi) ’ ( 3 ) 

n 

where V qn = V\ qn \, B = L/2it , the wavenumber q n is given 
by q n = nAq with A q = l/f2 and n = 0, ±1, ±2, 

We assume V qn is a continuous function at q n = 0 in the 
continuous limit A q —> 0 , and satisfies the condition, 

0(kn| 3/2 ) < I'M < 0(|<M| 1/2 ) (4) 

for q n —» 0 , in order to avoid a singular transport process 
characteristic in ID system (see (1671) for the decay rate, 
and (l 86 l) for the diffusion coefficient). 

In this paper, we consider the thermodynamic limit for 
the heavy particles, 

N 

L —> 00 , N 00 , c = — = finite, (5) 

1J 

where c is the concentration of the heavy particles. In this 
limit, we have A q —> 0 and the wavenumber and the mo¬ 
mentum become continuous variables. At an appropriate 
stage, we shall replace a summation with an integration 
and a Kronecker delta S Kr with a Dirac (5-function as 

A E ^ n ^S Kr (P - P') -> 6(P - P'), (6) 

1 q ' 

with Qn = tt/h (Hereafter we use a conventional notation 
'fPq for and drop the index n in q n ). 

Time evolution of the density matrix p(f) obeys the 
Liouvillc-von Neumann equation, 

= L np(t ), (7) 


where Lq is the unperturbed Liouvillian associated with 
Hq and gLy is the interaction Liouvillian associated with 
gV. In this paper, we focus our attention on the time evo¬ 
lution of the reduced density matrix for the test particle, 
which is defined by 

/(f) = Tr hev .[p(t)], (9) 

where Trhev. denotes the partial trace over the heavy par¬ 
ticles. 

We assume that the initial condition of the system is 
given in the form, 


P(0) = /(0) ® Phev.J (10) 

where p^f ev is the Maxwell distribution of the heavy par¬ 
ticles with temperature T, 


eq 

Phev. 


N 

n«p 

3=1 


Pi 


2 Mk B T 




Pi 


2 Mk B T 


(n) 

and ks is the Boltzmann constant. In the thermody¬ 
namic limit the time evolution of the density matrix as¬ 
sociated with the heavy particles is negligible since its 
deviation from is proportional to 1 /L. 

In order to discuss the “coordinate” and the “momen¬ 
tum” dependence of the distribution of the quantum par¬ 
ticles in parallel with classical mechanics, let us introduce 
the Wigner distribution function, which is a quantum 
analog of the phase space distribution function in classi¬ 
cal mechanics^, GS , defined by 


p w (X,{X j },P,{P j },t) 

V E PkAkM P APj},t)e i{kX+klXl+ ' +kNXN) , 


n N + 1 


M fc B 


( 12 ) 


with the Fourier component 


Pk,{kM P >{P3ht)=^/^ 


X 







(13) 


where the volume factor y in (fl3l) comes from the 

normalization in the representation of the Wigner basis 
(see (EfflO . Here each of the notation {A,}, {P 3 } repre¬ 
sents a set of variables for the N heavy particles, and 
the momentum states of N + 1 particles | p, {p.,}) is an 
eigenvector of Ho, 


H o\p, {Pj}) = £ P ,{ Pj }\P, {Pj}), (14) 


where Lh is the Liouvillian that is defined by the com¬ 
mutation relation with the Hamiltonian as 


LhP 



(Lq + gL v )p, 


with 


e p.{pjl — 


2 m 


N 

E 

3=1 


pj_ 
2 AT 


( 8 ) 


(15) 
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The “wavenumber” k and the “momentum” P in the 
Fourier component of the Wigner representation are 
given in terms of the momenta in the matrix elements 
{P, {Pj}\p(t)\p',{Pj}) as 

hk=p-p', P=)-{p + p'), 

2 j (16) 

hkj = Pj - p'j , P :j = - (pj + p'j). 

III. COMPLEX SPECTRAL REPRESENTATION 
OF THE LIOUVILLIAN 

A. Preliminaries 

In this subsection, we give a brief summary of the Liou- 
ville space representation of the Hilbert space operators, 
(see [4[ for more detail) 

The Liouville space is spanned by linear operators A, 
B } ■ ■ ■ in the ordinary wave function space. We represent 
these linear operators by double bra-ket notations as |A}}, 

|!?)},••• in order to distinguish them from the bra-ket 
notation in the wave function space. The inner product 
in the Liouville space is defined by 

«P|A» ee Tr [B + A], (17) 

where B + is the Hermitian conjugate of a linear operator 
B in the wave function space. 

In the Liouville space, we can consider a linear opera¬ 
tor P acting on vectors, which is operators in the wave 
function space. We call it a super-operator in case we 
want to emphasize the difference between it from the lin¬ 
ear operator in the wave function space. One can define 
a super-operator that is the Hermitian conjugate to 
T as 

((A\Ti\B)) = (({B\T\A)))°-°-, (18) 

where the notation c.c. means the complex conjugate. 
Then, one can introduce the Hermitian super-operator 
which satisfies 

T t = T, (19) 

and the unitary super-operator which satisfies 

t/l = U- 1 . (20) 

The Liouvillian is an example of the Hermitian super¬ 
operators in the Liouville space, i.e., 

L' h = L h , (21) 

and the time evolution operator 

U(t) = e~ iLat , (22) 

is an example of the unitary super-operators. 


We now introduce a basis set in the Liouville space for 
ID quantum Lorentz gas in terms of the eigenstates of 
Hq in JTU). Let us introduce an abbreviated notation 
of the N + 1 variables p = {p, pi,-- - , pn} to avoid 
heavy notations. Then, we express the eigenstate and the 
eigenvalue of H 0 in (fill) as |p) and e p . The eigenstates 
of Hq satisfy the orthonormality relation 

N 

<p|p'> = ^ r (p-p') = s Kr (p-p') n a* r c pj-p'j )> (23) 

i=i 

and the completeness relation 

EIpXpI = l ( 24 ) 

p 

We define vectors in the Liouville space by 

|p,p'» = IpXp'I- (25) 

These vectors are eigenstates of L 0 , 

A)|P,P')) = ^(e P - ep/)|p,p')). (26) 

The set of eigenstates {| p, p'))} forms a orthonormal ba¬ 
sis in the Liouville space (see (ED) 

«P, p'Ip"> P'")) = ^ r (P - P'O^^P' - P"0, (27) 

and the completeness relation 

EE Ip.p , )x<p.p'i = i - (28) 

p p' 

We note that the eigenvalues of L 0 in (1M1) are degen¬ 
erated (e.g. for p' = p) even though e p are not degener¬ 
ated. This is the intrinsic degeneracy of the Liouvillian 
mentioned in the introduction. 

In terms of the basis, we have a simple expression for 
a matrix element of an operator A in the wave function 
space as 

^p.p' = <P|-4|P'> = ((P, P'l^»- (29) 

A matrix element of the interaction Liouvillian is given 

by 

<<p,p'|^|p",p"')> 

= ^[gVp.p"5 Kr ( p' - p'") - 5 Kr {P - p")gV p i" t pif. ^ 

We define the Wigner basis by 

|k;P)) = ^+ r |p,p / )), (31) 

with the relations in CSD, where the volume factor 
fl n +1 is introduced to take a continuous limit —> oo 

on the variables P , {Pj} (see © and (15^1) 1 analogous 
to the classical system where momenta are intrinsically 
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continuous variables @. The set of the Wigner basis also 
forms a complete orthonormal basis 

«k;P|k , ;P')) = S Kr {k - k')Sn K {P - P') 

N 

(32) 

3 = 1 


with 


S nh (P-P') = n h S Kr (P-P'), (33) 

E?wEi k ; p Mk;P| = i. (34) 

k R p 


With the Fourier component of the interaction potential 
in 0, the Wigner representation of the matrix element 
of the interaction Liouvillian (13U1) is given by 


«k;P|( ? Lv|k / ;P , » 

1 N T/ JV_1 

= ^ YP^S Kr {k - k' + % - k') n S Kr (ki - K) 


.7=1 


x[7)| (fe k,) flp} (k fc,) 

xJnJP-P'), 


»(#j) 

„ ——fe') — fe')i 

Vp Vpl \ 


(35) 


with ftk' = p" —p'" and P' = (p" + p"')/2, where rf P and 
7 )p^ are displacement operators acting on the momenta P 
and Pj , respectively, as 


€>/0P) = /(P +p), V%f(Pj) = f(Pj +Pj )• (36) 


B. General formalism 

The eigenvalue problem of the Liouvillian is given by 

= ^IPM)), (37a) 

((pW|P if = ((pM|ZM, (37b) 

where the indices a and p specify an eigenvalue (espe¬ 
cially v denotes the spatial correlation subspace (see jj])), 
and | and ({F^\ are right- and left-eigenstates, re¬ 
spectively. 


We introduce effective Liouvillians by using the 
Brillouin-Wigner-Feshbach formalism 0,[l3 with projec¬ 
tion operators P^ and Q^ satisfying 

P^L 0 = L 0 P^\ 

(38) 

pWpbd = S ^, 

(39) 

=/*+!, 

(40) 


V 


and 


pM + QM=7 J v +1) (41) 

where In+i is the unit operator for the N + 1 particle 
system. By applying these projection operators on (I37al) , 
the eigenvalue equation of the Liouvillian takes the form 

^M(^H)pW| F M» = zMpM\FM)), (42) 


where 


¥"\z) = P {v) L h P M 

+ P^LuO^ _-_QOOp pO) 

+ Lh ^ z-QMLhQM^ H 

is the effective Liouvillian. Its second term is the self¬ 
frequency part that corresponds to the self-energy part 
of an effective Hamiltonian in the case of the Hamilto¬ 
nian operator in the wave function space. The effective 
Liouvillian is also called the collision operator which is 
a central object in the kinetic theory in non-equilibrium 
statistical mechanics HEP- One can see from its eigen¬ 
value equation (TT 21 ) that the collision operator shares the 
eigenvalues with the Liouvillian. The eigenvalue equa¬ 
tion of the collision operator (14211 is non-linear , i.e. the 
collision operator itself depends on the eigenvalue. 

In terms of the right- and left-eigenstates of the colli¬ 
sion operator (z), the right- and the left-eigenvectors 
of the Liouvillian L p are given by 

| fM» =[PM +CM(ZW)]PM|PM)), (44a) 
((pH| = ((pM|pM [pM + dM(zM)] , (44b) 

with the creation-of-correlation operator 

C {v) {z) =--4- ttQ 1v) LhP M , (45a) 

v ; z - QMLhQM 

and the destruction-of-correlation operator 

£>M(z) = pMLhQM -4- rT , (45b) 

v ’ z-QMLhQW 

which are off-diagonal transitions between the sub¬ 
space and the P^ subspace Q. 

It is well-known for an unstable quantum system with a 
continuous spectrum that the effective Hamiltonian be¬ 
comes a non-Hermitian operator due to the resonance 
singularity in the self-energy part |n|. Similarly, the col¬ 
lision operator becomes a non-Hermitian operator in the 
Liouvillc space in the thermodynamic limit. As a result, 
the collision operator has eigenstates with complex eigen¬ 
values that are called resonance states. The imaginary 
part of the complex eigenvalue of the Liouvillian gives 
a transport coefficient of the system [l3| (see also (I85al) 
and ». 
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C. Application to the weakly-coupled ID quantum 
perfect Lorentz gas 

Let us now apply the general formalism presented 
above to the weakly-coupled ID quantum perfect Lorentz 
gas. Using the Wigner basis m, the Fourier component 
of the Wigner distribution function (fldl) is represented 
by 

P k<{ki }{P,{Pj},t ) = ((fc,{%};P,{P,}|p(t))). (46) 

In order to apply the Brilloun-Wigner-Feshbach for¬ 
malism, we define the projection operators as 




(47) 


and 


Q(k) = 1 _p(k) = £ -_£|k ; P)X<k ; P|, (48) 


n N+1 

k^(fe,{0}) h P 

where we have used the notation {0,} to indicate that all 
wavenumbers associated to the heavy particles are zero. 
For the projection operator P^, we have 

l. p 

p( fc )LoP (fc) |/c,{0,};P)) = — |&,{0 j };P», (49) 


and 


gP^LvPW = 0, 


(50) 


because of Vo = 0 due to the assumption (HI). 

In the weak coupling situation, the collision operator 
can be approximated up to the second order in g as 


E \z) = P (fe) L 0 P (fe) 


+ g 2 P ( ' k) L v Q {k) -Q (fc) L y P (fc) . 


(51) 


z- L 0 

Using (1551) and (1551) . we have 

«M0,};P|44 fc) MIM0,};P'» 


N 




i= i 9 

1 


h 2 


z — (k — q)P/m — qPj/M 
x5n h ( P-P'). 


tiFfipr-tiptt:) 


(52) 


We focus our attention on the test particle. We define 
the reduced collision operator for the test particle as 


’EE) = T>hev.[^2 K EKevJ 


O'), 


(53) 


where we put a bar on the notation to distinguish it from 
the collision operator for the whole system (15TT) . Let us 


now derive an expression of its matrix element, 

«fc ; piEEp;n> 


= «*;P| 


L {pj} h 


\k;P'))- 

(54) 


Taking a large volume limit L —> oo for the momentum 
variables and performing the integration over the mo¬ 
menta with the aid of CD and m, we have 


«fc;P|EE)|fc;^'» = 


fcP _ 2 _j_ y- 
m ^ n 2 


dPj 


3= i' 


><E 

9/0 


in ' 2 


-^-('np 9 Vp^ q -Vp" q Vp^) 


z — (k — q)P/m — IPj /M 

XPhev.^lWP-P'), 


(Vp q Vp? q - Vp^ 9 Vp.) 


(55) 


where ( Pj ) is the Maxwell distribution for the heavy 
particle j with temperature T, 


PZM) 


1 




1/2 


27 rMkpT J 


exp 


P 2 

3 


2 Mk B T J 


V 


(56) 


The first and second terms in (1551) are called the “flow 
term” and the “collision term”, respectively. 

Since the mass of the test particle m is much smaller 
than the mass of a heavy particle M, we expand the 
propagator in (1551) as a power series of the ratio m/M. 
Then we have 

«*; ^lEE)!*; n = «fc; p’)) + o , 

(57) 

where 


((fe;P|#)(z)|fc;P')) = 


kP 2t xg 2 c 1 , 

m h 2 

9/0 


xd 


hq/2 


z — (k — q)P/i 


-a 


hq/2 


WP-P'), (58) 


is the collision operator for the perfect Lorentz gas 
(m/M -A 0). Here is defined by 


Z.hq/2 _ „§9 .-§9 

= Or ~ Vp ■ 


(59) 


Note that there is no temperature dependent term in 
(1551) . Temperature dependence comes from the correction 
term to (1551) starting with the first order in (m/M) (see 
Appendix). 

For the reduced collision operator, we write the eigen¬ 
value problem as 

E ) (E ) )lE ) » = 4 fe) lE ) », 


(60a) 
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((vM\^ k \zW) = Z M{(vW\. (60b) 

We note that = Z^ for our Lorentz gas, because the 
heavy particles are in an eigenstate with zero eigenvalue, 
i.e. they remain in thermal equilibrium. 

The collision operator in (l60l) depends on its eigen¬ 
value. In this sense, the eigenvalue equation is still non¬ 
linear. Our main goal is to construct the solution of 
the nonlinear eigenvalue problem. As a preparation to 
achieve this goal, however, we here restrict our atten¬ 
tion on the linear approximation. The complete basis 
obtained for the linear problem presented in this paper 
will be used to achieve the main goal in a forthcoming 
paper. 


IV. THE BOLTZMANN COLLISION 
OPERATOR 

To solve the eigenvalue problem (15U1) . we here study a 
situation where the wavenumber k satisfies 

\k\ < gT 1 , (61) 

where d is the interaction range between the particles. In 
this situation in addition to the weak-coupling, we shall 
show that the eigenvalue dependence of the collision op¬ 
erator (1581) is negligible and the collision operator if>( k \z) 
is reduced to the phenomenological Boltzmann collision 
operator. 

For a spatial inhonrogeneity satisfying the condition 
©, a typical value of q appearing in |3]) is much larger 
than k in (1551) . 


1*1 « M- (62) 

Then we can neglect k in the denominator in the second 
term in (1551) . We note that since the collision term in (1551) 
is an operator acting on the momentum P, this term 
does not commute with the flow term. This property 
changes the fc-dependence of the eigenvalue of from 
the linear ^-dependence in the flow term to some different 
fc-dependence. Indeed, the component of the eigenvalue 
associated to the diffusion processes starts with k 2 in the 
series expansion in k (see (185al) ). 

On the other hand, we may expect that the imaginary 
part of the eigenvalue z in (l60l) is proportional to g 2 for 
g <Cl because of the factor g 2 in front of the collision 
term in (1581) . If this is the case, we can replace z in (1551) 
by +*0, 


V> (fe) (4 fc) ) = ^ (fc) (+*«) + 0(g 4 ). (63) 

Here -MO means that the collision operator ^ k \z) is 
evaluated on the real axis approaching from the upper 
half-plane of z to ensure the time evolution is oriented 
to the future t > 0 Q. Combining these arguments, we 
can approximate ip^ k \z^' > ) by the new collision operator 


given by 


((fc;P|Vi fe) |fc;P')> = 

/ OO 

dq\V q \ 2 dp q/2 — 

-oo +* 

xS(P-P'), 


kp 

m 


2tt g 2 c 
1 


aM/2 


-Me + qP/r 


(64) 


where we have taken the thermodynamic limit d5j). This 
is the phenomenological Boltzmann collision operator for 
the ID quantum perfect Lorentz gas. M 

(k) 

We note that in spite of the fact that ip B is in the 
subspace that is orthogonal to the P (0 ^ subspace, 
the second term of the matrix element of ' does not 
depend on k and is the same as the matrix element of 
the i/j b ' > that is associated to the spatially homogeneous 
component k = 0 and generates a time evolution of the 
momentum distribution function. In this sense the phe¬ 
nomenological Boltzmann collision operator is a “coarse¬ 
grained object” in space, and it does not contain infor¬ 
mation of microscopic structure of the interaction. As a 
result, the Boltzmann collision operator is not applicable 
to a situation with so small scale of spatial structure that 
the microscopic shape of the interaction starts to play a 
role. 

We also note that the second term in (1551) vanishes for 
the classical ID perfect Lorentz gas. This is because the 
momentum distribution function cannot change in time 
by an elastic two-body collision in the classical ID sys¬ 
tem (for the perfect Lorentz gas, or for the system that 
consists of same mass particles). This is not the case in 
quantum system, because there is a forward scatting due 
to a quantum effect in addition to a backward scatter¬ 
ing. In this sense, the quantum system has an advantage 
on the mathematical simplicity over the classical corre¬ 
sponding systems. 

Furthermore, we note that the above derivation of the 
Boltzmann collision operator is still semi-quantitative, 
because the applicability of this approximation as m 
is too rough. Precise limitation of the Boltzmann ap¬ 
proximation will be clarified in the forthcoming paper by 
solving the nonlinear problem of the collision operator 
with the aid of the complete eigen-basis of the Boltzmann 
collision operator obtained in this paper. 

By performing the q integration in (1641) with a formula 


lim -= —in8(x) + V—, (65) 

e->+o x + ze x 


where V denotes the principal part, a matrix element of 
the collision operator is expressed as 

((k-,P\^ ) \k-,P , )) = —8(P-P') 
m 

2 

+i^^[5{P + P’) - 5{P ^ P% (66) 

with 


87 T 2 mc 
h 2 \P\ 



(67) 


7 p = 











We have yp —t 0 for P — > 0 (see (HD). 

From the expression (1551) . one can see that the Boltz¬ 
mann collision operator has non-vanishing matrix ele¬ 
ments only between the states | k; P)) and | k; — P)). Phys¬ 
ically, this is because there are only forward and back¬ 
ward scattering in the ID quantum system. Hence, in 
terms of this basis, the Boltzmann collision operator is 
represented by a 2 x 2 non-Hermitian matrix, 


/ (*0 

vv = 


kP/m — ig 2 "fp/2 ig 2 "fp/2 

ig 2 ^p/2 -kP/m - ig 2 7p/2 


( 68 ) 


(k) 

In terms of the operator %jj y B ', the time evolution equa¬ 
tion for the reduced density matrix for the test particle 
is approximated by 

i^ t P {k) \f{t))) = V4 fc) P (fc) l/W}}, (69) 

where 

p {k) = 7^Ei fc ; p »« fc ; p i- ( 7 °) 

“ft P 


V. EIGENSTATES OF THE BOLTZMANN 
COLLISION OPERATOR 

In this section, we present the eigenstates of the Boltz¬ 
mann collision operator (16811 and summarize their prop¬ 
erties. The structure of the collision operator (1551) is so 
simple that the solution of its eigenvalue problem itself 
has already been known M- We present here the explicit 
form of the solution in order to show a new physical in¬ 
sight into the time evolution generated by the Boltzmann 
collision operator, especially for the distribution function 
localized in moderately small spatial region of a scale less 
than the mean-free-length, but yet with a large enough 
width as compared with the microscopic scale given by 
the interaction range. In the following, we shall sepa¬ 
rately treat situations k = 0 and k ^ 0. 


A. The spatially homogeneous situation: k = 0 


8 

has a solution consisting of the eigenvalues and eigenvec¬ 
tors, 

4% = 0, (73a) 

with 

l ^ )) = 7f (|0;|jP|)) + |0;_|p|))) ’ (73b) 

and 

Z-:P = -ig 2 lP, (74a) 

with 

l0 -; )p))= 72 (|O;|p|))_|O;_ ™- (74b) 

Here we note that |</>±.p)) = |<£±?_p)) and z±) P = z±}_ p . 
The eigenstate |<^+.p)) is the equilibrium mode, while the 

eigenstate |0_.p)) is the decay mode of the momentum 
distribution function. As a result, the momentum relax¬ 
ation time Tp is given by 


1 

rp = ——. 

9~1p 


(75) 


Since the operator ED is a symmetric anti-Hermitian 
matrix, each left-eigenvector is the transpose (Hermitian 
conjugate) of the corresponding right-eigenvector (I73bl) 
and (I74bl) . 


«^pI = ^(«0; |P|| + ((0; -\P\\) = «E° ; U (76) 


u%\ = i=(«0; \P\\ - ((0; -\P\\) = U%\. (77) 

These eigenstates satisfy the bi-orthonormality, 

Ua-pl^-P')) = <W P(^ - P') + S(P + P')} , (78) 

and the bi-completeness relation, 

_ fOO 

E / (79) 

«=± 7o 


We first consider the spatially homogenous situation 
in the k = 0 subspace. This subspace corresponds to 
momentum distribution function. Then, the Boltzmann 
collision operator is given by an anti-Hermitian matrix 

< = (_\ Y) ■ (7D 

The representation ED is given with respect to the basis 
consisting of |0;P)) and |0; — P)). The right-eigenvalue 
equation 


= 4;pl^i°p)): 


where the integration over P is restricted to the region 
P > 0 to avoid double counting. 


B. The spatially inhomogeneous situation: k ^ 0 


Let us now consider the spatially inhomogeneous situa¬ 
tion with k ^ 0. In this situation, the eigenvalue equation 
of the collision operator (l68l) is 


det[v4 fc) 


zh 1 = 


■g 2 ip 


2 



g 2 ip ^j 


= o, 
(80) 


(72) 
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where I 2 is the unit matrix of size 2. Hence, the eigen¬ 
values are 


where 


(k) 

4 ;P 


■g 2 ip 


± 


\p\ 


m 


(k 2 -klf/\ 


g 2 ip = J_ 
2\P\/m Ip ’ 


(81) 

(82) 


is a wavenumber that is equal to the inverse of the mean- 
free-length Ip of the test particle with momentum P. 
Here we choose the branch of the square root function as 


f \Jk 2 — kp (fc > kp) 

(fc 2 - fcf ,) 1 / 2 = l iy/ k 2 P - fc 2 (|fc|<fc P ) (83) 

— \/fc 2 — kp (k <—kp) . 


These eigenvalues tend in the limit \k\ —> 0 as 


(fc) 

z +'p 

(fc) 

z -’p 


z (0) 

c +-,p-> 

z (°) 

— — ; P ? 


(84) 


with (173a!) and (174al) . 

For \k\/kp -C 1, one can expand the eigenvalue z^°} p 
and z^} p as 


(fc) ^ 

+;P “ 


z 


■ g 2 ip ( k \ 2 
1 4 \kp) 


and 


■g 2 ip 

16 



(85a) 



• 2 , - g 2 ip 

—ig 7 p + 1 —-— 


+0 



6 



.g 2 ipf k Y 
1 16 \k P ) 


(85b) 


First few terms of the expansion (185al) give transport 
coefficients of hydrodynamic equations. For instance, the 
first term, which is second order in k , gives the diffusion 
coefficient, 


= g 2 lP = (P/m) 2 

- 4 k 2 p g 2 jp 


( 86 ) 


We have Dp —> 0 for P —> 0 (see 0 ). The higher order 
terms of the expansion also give transport coefficients of 
the Burnett equation (see e.g. 0 )- 

In FIGH] we show a fc-dependence of the real part and 
the imaginary part of the eigenvalues for P 7 ^ 0. In the 
figures, the dashed lines and the dot-dashed lines rep¬ 
resent the eigenvalues z^) p and z^) p , respectively. The 
solid lines represent that these two lines overlap. Note 
that the eigenvelues (1511) can be rewritten as, 



g 2 ip 




1/2 


- 1 


(87) 



FIG. 1. Eigenvalues of the Boltzmann collision operator m 
are drawn as functions of k. (A) is the real part and (B) 
is imaginary part. In each figure, the dashed lines represent 
eigenvalue with a = + and the dot-dashed lines represent 
eigenvalue with a = —. The solid lines represent that these 
two lines are overlapping. The gray lines in (A) represents 
eigenvalues of the Liouvillian for a free light-mass particle 

y = ±(i/ 2 )(fc/fcp). 


Hence the eigenvalues depend on P only through the P- 
dependences of g 2r )p and kp. As a result, the spectrum 
is universal in the sense that the eigenvalues are inde¬ 
pendent of P, if we use the P-dependent units in which 
g 2 jp = 1 and kp = 1. 

The eigenstates corresponding to the eigenvalues in 
(18T1) are 


1 *£>,» = 


1± 


(t 2 -4)7 21 ' /2 


.-i*i 


IT 


k 

(fc 2 — fcp) 1 / 2 


I fc; \p\)) 

1/2 


( 88 a) 


\k;-\P\h 


((£ 


gpi 


\_k\_ 

k 


i± 


(fc 2 — fcp) 1 / 2 


1/2 


« fc ; 1^1 


T i 


1 T (fc 2 ~ k 2 P ) 1/2 
k 


1/2 


( 88 b) 


«fc;-|P| 


respectively. We note that z±) P = z±/_ p and |x± p)) = 
|x±;-p))- Moreover, we note that ((x± pl is n °t the Her- 


mitian conjugate of |x±.p}), because i/p’ with fc 7 ^ 0 is 


(fc) 
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neither a Hermitian operator nor an anti-Hermitian op¬ 
erator. 

In (|88l) . we have not yet normalized the eigenvectors 
by taking account of the fact that we have a diverging 
normalization constant at k = ±kp (see (1911) 1. The inner 
products of these right- and left-eigenstates are given by 


<«gpixS»>> = ± 


2(k 2 - kl) 1 / 2 

W\ 


[6(P - P') + S(P + P')}. 

(89) 


Then, normalized eigenstates for k ^ ±kp are given by 

i4‘W> = \Agtixg,)), 


sg’pi e «*£>,i, 

where the normalization constants are 

\k\ 


N ±-p = ± 


2(k 2 ~k 2 P y/ 2 ' 


(90) 


(91) 


For k ^ Pkp, they satisfy following bi-orthonormal and 
bi-completeness relations, 


the usual degeneracy point of the eigenvalues of a Hermi¬ 
tian operator for which a degenerate eigenvalue is shared 
by two distinct eigenstates. Such a branch point of the 
spectrum of a non-Hermitian operator in the parameter 
space are called an exceptional point (EP) [16}, which is 
also known as a non-Hermitian degeneracy points [13 . In 
the previous paper 0 , we have introduced a divergence 
free representation at EPs by continuously extending the 
Jordan block representation away from EPs. 

(k) 

Note that the imaginary part of the eigenvalues z±. P 
remain constant —g 2r )p/2 for |fc| > kp , and they do not 
vanish in the limit |fc| —> oo. As a result, one cannot 
get any information of the microscopic space structure of 
the interaction through the Boltzmann collision operator. 
This is clearly a limitation of the Boltzmann approxi¬ 
mation to analyze microscopic origin of irreversibility in 
dynamical processes. 


VI. TIME EVOLUTION OF THE WIGNER 
DISTRIBUTION FUNCTION AND THE 
EXCEPTIONAL POINT 


«#$>!<^P'» = <*«,«' [S(P - P') + S(P + P')} , (92) 

/>oo 

£/ dP\<f>%n$$p\=P {k) - ( 93 ) 

a =± J ° 


The points k = ±kp in Ea. flSU) are branch points of 
the eigenvalues in the complex k plane. At these points, 
both eigenvalues and eigenvectors coalesce as, for k = kp, 


(k P ) .g 2 ip 

Z ±:P = 1 2 ’ 

(94a) 

|x±.p » = \kp; |P|)} + i\kp-, -|P|)), 

(94b) 

{(x ( ±^\ = ((kp-,\P\\+i{(k P -,-\P\\, 

(94c) 

and as, for k = —kp, 


^ 3S- 

II 

1 

<s>. 

to ^ 

"a 

(95a) 

|x4“p p) » = | - fc P ; |P|)) - - k P ; -|P|», 

(95b) 

«3cfcp p) l = -«-*p; 1^11+ i«-*p;-1^11- 

(95c) 

Since there is only one linearly independent eigenvector 
at the points, the Boltzmann collision operator (l68l) can 
not be diagonalized. Instead, the collision operator has 
the Jordan block structure at these points (see fl5l Il8}). 
This coalescence of eigenvectors does not take place at 


In this section, we discuss the relation of the spectral 
structure to the transport processes in the system. In the 
previous paper 0 , we have shown that the exceptional 
points in the spectrum of the Boltzmann collision opera¬ 
tor lead to the telegraph equation in the time evolution 
in coordinate space. Here, we study further detail of the 
time evolution. The results presented in this section are 
summarized as follows: 

We found shifting motion of the peak of the distribution 
function in space in addition to spreading as a result of 
a diffusion type process. However, the mechanism of the 
shifting motion is very different in the domain \k\ < kp 
from that in the domain |fc| > kp. 

1. For |fc| < kp, the shifting motion comes from asym¬ 
metry in the momentum distribution before the 
momentum relaxation is complete. 

2. For |fc| > kp, the shifting motion comes from the 
real part of the eigenvalue that leads to a wave 
propagation with the initial velocity P/m. 

To see this, we shall analyze time evolution of the 
Wigner distribution function for the test particle 

/ OO 

dk f k (P,t)e ikX , (96) 

-OO 

with 

f h (P,t) = ((k;P\m)). (97) 

Multiplying ((k; P\ and using the completeness relation 
(1U51) to the formal solution of the Boltzmann equation 

P w \m)) = e-*™ t pW\m)), 


(98) 
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we have the following expression of the Fouier transfor¬ 
mation of the Wigner distribution function fk(P,t) 


fk(P,t) = 



The first two terms correspond to time evolution due 
to forward scattering and the third term corresponds to 
time evolution due to backward scattering. 

Note that the expression in (l99l) is regular at the EPs 
where z+) p = z^} p . Indeed, by taking the limit k —> kp 
or k —> —kp, we have 


f±k P (P,t)=e~ L ^ t f ±kp (P, 0) 

k P P . ( 100 ) 


a~ip + 

=F i - te 


m 


IJZte-^UU-P, 0 ), 


which is identical with the expression derived in terms of 
the Jordan block representation presented in 0. 

In the second and the third terms in (11001) . we have 
critical damping factors t exp[— (<? 2 7p/2)t]. Since the 
eigenvalues z+) P and z^) P change their values from pure 
imaginary values for |fc| < kp to complex values for 
|fc| > kp, the time evolution of the Fourier component 
(l99]l changes from an over-damping to a damped oscilla¬ 
tion as a critical damping behavior at the EPs. 

In order to see the relation of the spectral structure to 
the time evolution of the Wigner distribution function, 
we consider the following situation as an initial condition, 



- t = 0 . t = 1/(5 2 7p) 

. t = 5/ {g 2 1p) . t = 10/(5 2 7p) 

FIG. 2. Time evolution of the Wigner distribution function 
for P > 0. The solid line represents the initial distribution. 
The initial distribution is given by Eq. USD with kb = kp. 
The initial distribution evolves to the distributions repre¬ 
sented by the dotted lines as represented by the arrows. This 
figure does not depend on the value of P. 


We use the units in which Ip = 1 /kp = 1 and 
Tp = 1 /(<? 2 7 p) = 1, when we present results of numeri¬ 
cal calculations (see e.g. FIG®. With these units, the 
eigenvalues and the eigenvectors are independent of the 
value of P (see (1571) and (1551) '). Therefore the dynamics 
described by the Boltzmann equation is universal in the 
sense that the time evolution equation in the units is the 
same regardless of the actual value of P. 


A. Time evolution with the spectrum in |fc| < kp 

Let us first consider the situation where the initial dis¬ 
tribution is composed of the Fourier components with k 
in the region |fc| < kp. Hence, we take 


fk(P,0)= X k b (k)h(P), (101) 


kb < kp. 


(105) 


where Xk b {k) is a step function which is defined with a 
given value of kb by 


In this case, the time evolution of the Wigner distribution 
function for P > 0 is expressed by 


Xk b (k) 


1 (|fc| < k b ), 
0 (|fc| > k b ), 


( 102 ) 


and h(P) is a momentum distribution function that is 
normalized as 



dP h(P) = 1. 


(103) 


To extract the essence of the mechanism of the shifting 
motion, we here assume 


h(P) = 0 for P < 0, (104) 


i.e., initial distribution is composed of particles with pos¬ 
itive momentum. 


f w (X,P,t)= [^ dkf k (P,t)e ikX 

J-k h 


where 


k b 

kb 


dk- 


kp 


-k b 2 sjk 2 p - k 2 

e _l2 + ; p* cos (kX — <fk,p) 


■ O) 

+ e iz -,p 1 cos(kX + ipk,p) 


(106) 


fk(P, 0), 


Pk,p = arctan 


y/lej, - fc 2 


(107) 
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In FIGH we present the time evolution of Eg. (11 061) in 
X space. In the figure, the solid line is the initial distri¬ 
bution, and the dashed lines are the distribution of later 
times. Here, and in the following figures in this section, 
we have not specified the value of P, since the time evolu¬ 
tions presented in this figure and in the following figures 
are independent of P in the units in which Ip = 1 and 
t p = 1 . 

As shown in the figure, the peak of the distribution 
function shifts toward A' = Ip in the first stage of its 
time evolution. Afterwards, it stops the shifting motion 
and spreads its width symmetrically around the position 
X = Ip, where the time evolution is described by the 
diffusion equation with the diffusion coefficient Dp in 
(l86l) . This can be seen by the fact that (11061) can be 
approximated after the above mentioned first stage as, 

/ OO i 

dk-e~ Dpk2t cos [k(X - l P )]f k (P, 0) 

-CO " 

(108) 

which is a solution of the diffusion equation with a diffu¬ 
sion coefficient D P given by (15(1 . 

We note that, for t —> oo, the momentum distribution 
function 

1 Z 100 

f 0 (P,t) = —j dXf w (X, P, t), (109) 

is stationary as 

fo(P,t) = M-P,t) = i/o(P,0). (110) 


(A) t = 0 



This implies that the momentum relaxation has been 
completed, and this is the reason that the peak of the 
distribution no longer moves. 


Mechanism of the shift for the component with |fc| < kp 


In the above, we have shown that the peak of the distri¬ 
bution function in A' space shifts in spite of the fact that 
there is no real component of the spectrum for |fc| < kp. 
We now show that this is a result of the fact that the com- 
ponents of the distribution function correspond to z+-p 0 

and z^) P branches of the eigenvalues (I5T1) have asymme¬ 
try in X space, and their decay rates are different as 

|Im[z+.p]| < |Im[zi fe ; p]| for |fc| < kp. (Ill) 

(see also FIG.1(B)) 

To see this, let us first define the following two projec¬ 
tion operators 


poo 

= / dP\^ P m%\. 

Jo 

(112) 

*(k) ^(fe) i ^(fe) 

P K J = V+ +P- A 

(113) 


- f w (X, P,t) 

. f?(X,P,t) 

- fW P. t) 

FIG. 3. Time evolution of the distribution function 
f n ( X,P,t ), its “+” component fA and its ” component 
are presented for (A) t = 0, (B) t = 1/(<? 2 7p) and (C) 
t = 10/(g 2 jp). In each figure, the solid line represents / , 

the dashed line represents /^ and the dot-dashed line rep¬ 
resents /_ . The initial condition is given by Eq. ifToTTi with 
k = kp. This figure does not depend on the value of P. 


as in (1551) . 

Time evolution of both p y + ' and p_' components of 
| f(t))) with an initial condition (11(111) and with kb < kp 
for P > 0 are given by 


f± (X,P,t) = I dk((k-P\p^\f(t)))e 

J — OO 

/ k b 

dk __^ 

-k b 2y/tip — k 2 

(fc) 


ikX 


kp 


(114) 


x e lz ±-,p t cos(fcA =F <Pk,p)fk(P, 0). 


They satisfy 
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Then, we have 

f w (x, P, t) = /f (X, P, t) + fY (X, P, t). (115) 

In FIG.3, we present the time evolution of /T, fY and 
/" in X space with a fixed value of momentum P > 0. 
In each figure in the FIG.3, the solid line represents f w , 
the dashed line represents /Y and the dot-dashed line 
represents fY . As shown in the figures, the position of 
the peak of /Y is in the region X > 0 and the position of 
the peak of fY is in X < 0. As time passes, the peak of 
fY decays more rapidly than the peak of f'Y , and thus 
the peak of the total distribution f w shifts forward as 
shown in FIGS.3(B) and 3(C). 

Physically, the shifting motion corresponds to center 
of mass motion of the distribution due to asymmetry of 
the initial momentum distribution function. 


B. Time evolution with the spectrum in |fc| > kp 


In Sec. V, we showed that the spectrum of the Liouvil- 
lian takes complex values in the region |fc| > kp and their 
real part approaches to the eigenvalue for a free particle 
in the limit |fc| —> oo (see FIG.l). In this subsection, we 
discuss how this spectrum affects the time evolution of 
the system. For this purpose, we analyze time evolution 
with an initial condition (11011) with 


kb > kp. (116) 

For this initial condition, the time evolution of the 
Wigner distribution function for P > 0 can be divided 
into two parts; (i) fY with pure imaginary eigenvalues 
and (ii) f'Y with complex eigenvalues as 


/ OO 

dkf k (P, t)e ikx 

-OO 

= /f(X,P,t) + /f(X,P,f), 


(117) 


where (i) 


ff (X,P,t) = [^ dkfk(P,t)e 

J — kp 

kp 


ikX 


L 


k P 

kp 


dk- 


k P \J kp — k 2 

. O) * 

e lz +,p t cos(fcX — <pk,p) 

■ ( fe ) t 

+ e cos(fcX + ip k , P ) 


(118) 


fk(P, 0), 




Total 

\k\ < kp 
\k\ > kp 


FIG. 4. Time evolution of the Wigner distribution function 
for P > 0. The initial distribution (A) is given by (I101D with 
kb = 5.0 kp. The distribution functions are shown for (A) 
t = 0, (B) t = 1/(<? 2 7p) and (C) t = 5/(g 2, yp). In each figure, 
the solid line represents total value of the function f w , the 
dashed line and the dot-dashed line represents fY and fY , 
respectively. This figure does not depend on the value of P. 


and (ii) 


fY(X,P,t)=( [ " + f b ) dkfk (P, t)t 

\J—kb J kp J 

x e ik(,X-^Y 1 ~( k P/ k ) 2t ) _|_ 


AkX 


\J 1 - ( k P /k) 2 J 

i e *fe(A’+£y' 1 -(fcp/fc) 2 t) Y _ 1 _ 

V Y 1 - (fcp/fc ) 2 J. 

x f k (P, o). 


VH9) 
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Eg. (11181) consists of the over-damping components of 
m, while (11191) consists of the damped oscillation com¬ 
ponents. 

In FIG.4, we show the time evolution of and as 
well as the total distribution function f w with a specific 
value of momentum P > 0. In each figure, the solid 
line represents the total distribution f w , the dashed line 
represents fY and the dot-dashed line represents . 

The time evolution of fY has the same character as the 
time evolution of (11061) discussed in the above subsection 
A, namely the distribution fY ( X , P, t) firstly shifts to¬ 
ward X = Ip, within the time interval 0 < t < 1/(<? 2 7p). 
Afterwards, the distribution no longer shift and spreads 
its width by the diffusion process with its center fixed 
at Ip. On the other hand, the distribution f^(X,P,t) 
propagates as a wave packet with a velocity nearly equal 
to the initially give velocity P/m and decays in time. The 
wave propagation is due to the real part of the eigenval¬ 
ues of the Liouvillian in the region |fc| > kp. 


of the distribution function that spreads in moderately 
small spatial scale where hydrodynamic approximation is 
not applicable. There, we found two completely different 
mechanisms of the shifting motion; one is due to different 
value of the imaginary part of the eigenvalue of the Liou- 
villian, and dominates in a moderately small wavenumber 
associated to the inhomogeneity, while the other is due 
to the wave propagation associated to the real part in 
the eigenvalue that appears only in a moderately large 
wavenumber. 

In this paper we have not yet presented the solution of 
the nonlinear equation for the eigenvalue problem of the 
collision operator. Only by solving the nonlinear problem 
we can describe the irreversible processes occurring in mi¬ 
croscopic scale. The precise limitation of the Boltzmann 
approximation is also analyzable only after we solve the 
nonlinear problem. We will show the solution of the non¬ 
linear problem in a forthcoming paper where we will use 
the completeness relation of the eigen-basis of the Boltz¬ 
mann equation that is presented in this paper. 


VII. CONCLUDING REMARKS 

We have derived an irreversible kinetic equation for the 
weakly-coupled ID quantum perfect Lorentz gas through 
the spectral analysis of the Liouvillian starting from 
the first principle of physics without relying upon phe¬ 
nomenological arguments. As demonstrated in this pa¬ 
per, irreversibility is an intrinsic dynamical property of 
the basic laws of physics coming from the resonance phe¬ 
nomena that is associated to a complex eigenvalue of 
the Liouvillian (which corresponds to a pole of the re¬ 
solvent operator for the Liouvillian in the complex fre¬ 
quency space). In other words, the irreversibility is not 
an approximate concept coming from the human’s lim¬ 
itations on controlling microscopic information of large 
systems as often stated in many textbooks in statistical 
mechanics (see e.g. 0 )- 

Then, we have discussed a limitation of the applica¬ 
bility of the phenomenological Boltzmann equation. As 
shown in this paper, one cannot analyze the irreversible 
processes by the Boltzmann equation in so small a scale 
where the microscopic structure of the interaction be¬ 
tween the particles starts to play a role, because detailed 
space structure of the interaction is smeared out in the 
Boltzmann approximation. In order to describe the irre¬ 
versible processes in such a small scale, we have to solve 
the nonlinear eigenvalue problem for the collision opera¬ 
tor (1581) beyond the Boltzmann approximation. 

Moreover, we have shown that the perfect Lorentz gas 
with m/M —> 0 cannot deal with the temperature de¬ 
pendence. Hence, it is worth remarking that, in this 
approximation, one cannot study temperature depen¬ 
dent phenomena such as a transport process in the non- 
equilibrium steady state. In order to deal with the tem¬ 
perature dependence, we need to go beyond the lowest 
order approximation in m/M expansion. 

We have also discussed an interesting shifting motion 
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Appendix: Temperature dependence 


In this appendix, we give an expression for the cor¬ 
rection to the collision operator of the perfect Lorentz 
gas. 

We expand the propagator in (l55l) as a power series of 
the ratio m/M , 


1 

z ~ (k — l)P/m — IPj/M 

Y" (")" , A n 

^[z-(k-l)P/m] n + l \MJ ' 1 ’’ 

For each order term, one can perform the integration 
over momentum of N heavy particles using the following 
formulae for the Gaussian integrals, 
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where dPf 1 stands for the integration over momenta 

of N heavy particles, and 
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which is a result of the normalization. Note that the T 
dependence of (1551) comes from even power in the series 
expansion with ( m/M) in (IA.1I) . Hence, we 
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for the integrals in (1551) over the momenta of the heavy 
particles. As a result, T dependence appears already in 
the first order correction in m/M in spite that (l55l) has 
a contribution in its even power. Then, we have (1581) . 
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